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Abstract 

It appears that the development of machine vision may benefit from a detailed 
understanding of the imaging process. The reflectance map, showing scene radiance as a 
function of surface gradient, has proved to be helpful in this endeavor. The reflectance 
map depends both on the nature of the surface layers of the objects being imaged and the 
distribution of light sources. Recently, a unified approach to the specification of surface 
reflectance in terms of both incident and reflected beam geometry has been proposed. 
The reflecting properties of a surface are specified in terms of the bidirectional 
reflectance-distribution function (BRDF). 

Here we derive the reflectance map in terms of the BRDF and the distribution of 
source radiance. A number of special cases of practical importance are developed, in 
detail. The significance of this approach to the understanding of image formation is 
briefly indicated. 
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§1. THE REFLECTANCE MAP 

The apparent “brightness” of a surface patch depends on the orientation of the 
patch relative to the viewer and the light sources. Different surface elements of a non- 
planar object will reflect different amounts of light towards an observer as a 
consequence of their differing attitude in space. A smooth opaque object will thus give 
rise to a shaded image, one in which brightness varies spatially, even though the object 
may be illuminated evenly and covered by a uniform surface layer. This shading 
provides important information about the object’s shape and has been exploited in 
machine vision [1-8]. 

A convenient representation for the relevant information is the “reflectance map” [4, 
6]. The reflectance map, R(p, q), gives scene radiance as a function of surface gradient 
(p , q) in a viewer-centered coordinate system. If z is the elevation of the surface above 
a reference plane lying perpendicular to the optical axis of the imaging system, and if x 
and y are distances in this plane measured parallel to orthogonal coordinate axes in the 
image, then p and q are the first partial derivatives of z with respect to x and y. 

p — dz/dx and q = dz/dy 
* . 

The reflectance map is usually depicted as a series of contours of constant scene 
radiance (Fig. 1). It can be measured directly using a goniometer-mounted sample, or 
indirectly from the image of an object of known shape. Alternatively, a reflectance 
map may be calculated if properties of the surface material and the distribution of light 
sources are given. One purpose of this paper is to provide a systematic approach to this 
latter endeavor. Another is to derive the relationship between scene radiance and image 
irradiance in an imaging system. This is relevant to machine vision since gray-levels are 
quantized measurements of image irradiance. 

§2. MICROSTRUCTURE OF SURFACES 

When a ray of light strikes the surface of an object it may be absorbed, transmitted, 
or reflected. If the surface is flat and the underlying material homogeneous, the 
reflected ray will lie in the plane formed by the incident ray and the surface normal 
and will make an angle with the local normal equal to the angle between the incident 
ray and the local normal. This is referred to as “specular,” “metallic,” or “dielectric” 
reflection. Objects with surfaces of this kind form virtual images of surrounding 
objects. 

Many surfaces are not perfectly flat on a microscopic scale and thus scatter parallel 
incident rays into a variety of directions (Fig. 2a). If deviations of the local surface 







FIGURE 1 : A typical reflectance map for a surface, with both a glossy and a matte 

component of reflection, illuminated by a point source. The coordinates are surface 

slope in the x and y directions, and the curves shown are contours of constant scene 
radiance* 
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normals from the average are small, most of the rays will lie near the direction for ideal 
specular reflection and contribute to a surface “shine” or “gloss.” 

Other surface layers are not homogeneous on a microscopic scale, and thus scatter 
light rays which penetrate the surface by refraction and reflection at boundaries 
between regions with differing refractive indices (Fig. 2b). Scattered rays may re- 
emerge near the point of entry with a variety of directions and so contribute to 
“diffuse,” “flat,” or “matte” reflection. Snow and layers of white paint are examples of 
surfaces with this kind of behavior. Frequently both effects occur in surface layers, 
with some rays reflected at the nearly flat outer surface of the object, while others 
penetrate deeper and re-emerge after multiple refractions and reflections in the inhomo¬ 
geneous interior. 

The distribution of reflected light in each case above depends on the direction of 
incident rays and the details of the microstructure of the surface layer. Naturally, what 
constitutes microstructure depends on one’s point of view. Surface structures not 
resolved in a particular imaging situation are taken here to be microstructure. When 
viewing the moon through a telescope, for example, smaller “hillocks” and “craterlets” 
are part of this microstructure. This consideration leads to more complicated models of 
interaction of light with surfaces than those discussed so far. It is possible, for instance, 
to consider an undulating surface covered with a material which in itself already has 
complicated reflecting behavior (Fig. 2c). 

Reflectance is not altered by rotating a surface patch about its normal when there is 
no asymmetry or preferred direction to either the pattern of surface undulations or the 
distribution of sub-surface inhomogeneities. Many surface layers behave this way and 
permit a certain degree of simplification of the analysis. Exceptions are such things as 
diffraction gratings, iridescent plumage and the mineral called “tiger eye.” These all 
have a distinct directionality in their surface microstructure and will not be considered 
here any further. 

Considerable attention has been paid to the reflective properties of various surface 
layers. Some researchers have concentrated on the experimental determination of 
surface reflectance properties [9-21]. At the same time, many models have been 
developed for surface layers based on some of the considerations presented above 
[22-35]. Models often are too simple to be realistic, or too complicated to yield 
solutions in closed form. In the latter case, Monte-Carlo methods can be helpful, 
although they lead only to numerical specification of the reflecting behavior. Purely 
phenomenological models of reflectance have found favor in the computer graphics 
community [36, 37, 38]. Several books have appeared describing the uses of reflectance 
measurements in determining basic optical properties of the materials involved [39, 40, 
41]. Attention has been paid, too, to the problem of making precise the definitions of 

reflectance and related concepts [42, 43]. 
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FIGURE 2c: Compound surface illustrating more complex 
model of interaction of light rays with surface microstructure. 
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§3. RADIOMETRY 

# 

A modern, precise nomenclature for radiometric terms has been promoted by a 
recent NBS publication [43]. The following short table gives the terms, preferred 
symbols and unit dimensions of the radiometric concepts we will have occasion to u se 
for the development presented here: 


RADIOMETRIC CONCEPTS 

[W] 

[W.sr-1] 

[W. m-2] 

[W. m-2] 

[W.m-2.sr-l] 

<* 

Radiant flux, 4>, is the power propagated as optical electromagnetic radiation and is 
measured in watts [W], The radiant intensity, /, of a source is the exitant flux per unit 
solid angle and is measured in watts per steradian [W • sr -1 ]. The total flux emitted by 
a source is the integral of radiant intensity over the full sphere of possible directions 
(47r steradians). The irradiance, E, is the incident flux density, while radiant exitance, 
A/, is the exitant flux density, both measured in watts per square meter of surface 
[W • m -2 ]. The total radiant exitance equals the total irradiance if the surface reflects 
all incident light, absorbing and transmitting none. 

The radiance, L , is the flux emitted per unit foreshortened surface area per unit solid 
angle. Radiance is measured in watts per square meter per steradian [W • m -2 • sr -1 ]. It 
can equivalently be defined as the flux emitted per unit surface area per unit projected 
solid angle. Radiance is an important concept since the apparent brightness of a surface 
patch is related to its radiance. Specifically, image irradiance will be shown to be 
proportional to scene radiance. 

Radiance is a directional quantity. If the angle between the surface normal and the 
direction of exitant radiation is 0, then the term “foreshortened area" stands for the 

actual surface area times the cosine of this angle 0. Similarly the “projected solid angle" 
stands for the actual solid angle times the cosine of the angle 0. Here we will use the 
symbol w to denote a solid angle, while will be used to denote a projected solid angle. 


Radiant flux 




Radiant Intensity / = d$/du> 


Irradiance 


E = d*/dA 


Radiant Exitance M = d$/dA 


Radiance 


L = d 2 $ / (dA • cos 0 • did) 
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corresponding infinitesimal solid angles and projected solid angles 


dSl = du> • cos 6 

The following example (Fig. 3) will illustrate some of these ideas. Consider a source 

of radiation with intensity I in the direction of a surface patch of area dA, oriented 

with its surface normal making angle 6 with the line connecting the patch to the source. 

In fact, as seen from the source, it appears only as large as a patch of area dA • cos 0 

would when oriented perpendicular to this line. The corresponding solid angle is simply 

the area of this equivalent patch divided by the square of the distance from the source 
to the patch. Thus, 
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It dw and dQ are 
respectively, then 


du) es dA * cos 0 / r 2 [sr] 


The flux intercepted then is 


</#=/• du> *= / . dA • cos 6 / r 1 [W] 

The irradiance of the surface is just the incident flux divided by the area of the surface 
patch. 

E = d$/dA — I ♦ cos 6 / r 2 [W. nr 2 ] 


SURFACE 

NORMAL 





FIGURE 3 : Point source illuminating a surface, 
illustrating basic radiometric concepts. 
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§4. THE BIDIRECTIONAL REFLECTANCE-DISTRIBUTION FUNCTION 

The Bidirectional Reflectance-Distribution Function (BRDF) was recently introduced 
by Nicodemus, Richmond, Hsia, Ginsberg, and Limperis [43] as a unified notation for 
the specification of reflectance in terms of both incident- and reflected-beam geometry. 
The BRDF is denoted by the symbol f r and captures the information about how 
“bright” a surface will appear, viewed from a given direction, when it is illuminated 
from another given direction. To be more precise, it is the ratio of reflected radiance 
dL r in the direction towards the viewer to the irradiance dEj in the direction towards 
a portion of the source. In symbols, 



Here, 6 and 0 together indicate a direction, the subscript / denoting quantities 
associated with incident radiant flux, while the subscript r indicates quantities associated 
with reflected radiant flux [43]. 


The geometry is as depicted in the figure (Fig. 4). A surface-specific coordinate 
system is erected with one axis along the local normal to the surface and another 
defining an arbitrary reference direction in the local tangent plane. Directions are 
specified by polar angle $ (colatitude) measured from the local normal and azimuth 
angle 0 (longitude) measured anti-clockwise from the reference direction in the surface. 
In general, incident flux may arrive from many portions of extended sources, so incident 
radiance L/(0/,0/) is a function of direction. If we consider the component of flux </$/ 
arriving on the surface patch of area dA from an infinitesimal solid angle duij in the 
direction (0,-,0 z ) we obtain 

d$j ss Lj • cos 6j • dm} • dA = Lj • dQj * dA — dEj * dA 

where dEj=LfCos6j»dwj is the incident irradiance contributed by the portion of the 
source found in the solid angle du)j in the direction (0/,0/). Similarly, the radiant flux 
emitted into an infinitesimal solid angle d u> r in the direction (0 r ,0 r ): 


<P$ r — dL r ♦ cos 6 r • du> r • dA — dL r • dQ. r • dA 


where dL r (6 r ,0 r ) is the radiance in the direction (6 r ,<p r ) due to the reflection of the 
flux incident from direction (0 f -,0 z -). The notation dX, where X is one of the 
radiometric quantities introduced in §3, will always denote a directional quantity, that is, 
one which depends on either the incident or exitant direction. The notation cP-X will 
mean a bidirectional quantity, which depends on both the incident and exitant 
directions. Thus, the incident flux d$j depends only on the direction of incidence, but 

the exitant flux cP$ r depends on both the direction of emittance and (implicitly) on the 
direction of incidence. 
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<Pr 


FIGURE 4: Local geometry of incident and 
reflected rays needed for the definition of the 
bidirectional reflectance-distribution function 
(BRDF). Redrawn from [43]. 



From these values of incident and exitant flux, the BRDF is defined as follows: 

*r-0r) = (d** r /dQ r ) / </*,- = dL r /dE t 

and thus has dimension inverse steradian [sH]. The BRDF allows one to obtain 

reflectance for any defined incident and reflected ray geometry simply by integrating 
over the specified solid angles [43]. 


§5. INTEGRALS OVER SOLID ANGLES AND PROJECTED SOLID ANGLES 

The admitting aperture of an imaging system may occupy a significant solid angle 
when seen from the point of view of the objects being imaged. We will furthermore 
have to deal with extended sources. In both cases it is necessary to integrate various 
quantities over solid angles or projected solid angles. This can be accomplished by 




Horn & Sjoberg 


October, 1978 


double-integration with respect to the polar and azimuth angles (Fig. 5). 
quantity to be integrated, we have 

J w X dw = J I X sin 0 d0 d<t> 


If X is the 


and 


J r fn rn/2 

q X dQ. = J ^ J X cos 6 sin 0 d0 d<f> 


If for example X 
surface, then 


1 and the region of integration is the hemisphere above the object’s 

J * r * /*"/2 

dm — I I sin 0 rftf = 2>r 


while 


J r pir pn/2 

Q X dQ = J J 0 (1/2) sin 20 d0 d<f> = tr 


The latter result will be used in the discussion of perfectly diffuse reflectance. 
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§6. PERFECTLY DIFFUSE REFLECTANCE 

A perfectly diffuse or “lambertian” surface appears equally bright from all directions, 
regardless of how it is irradiated, and reflects all incident light [43]. Thus the reflected 
radiance is isotropic, that is, L r is constant, with the same value for all directions 
(d r ,0 r ). Also, the integral of reflected radiance over the hemisphere above the surface 
must equal the irradiance E. This implies that the BRDF for this ideal surface, f r id , is 
constant, and that the radiant exitance M equals the irradiance E. If the reflected 
radiance is L r , then the radiant exitance can be found by integration: 

M = Jq L r dSl r sss L r w 

From this one finds that 

f.r,id = L r /Ei = L r /M = \/ir 

If we have an extended source with radiance Lj , then the irradiance on the surface due 
to a small portion of solid angle dojj lying in the direction (0/,0/) is 

dEj ( 6j, 0,) = Lj(6j, 0/) cos 6j duij. 

The reflected radiance is then 

L r = (l/*r) cos 6j du>j 

This is a form of Lambert’s cosine law. 


§7. COLLIMATED SOURCES AND THE DIRAC DELTA-FUNCTION 

» 

Not all sources are extended. One way to deal with sources that are highly 
collimated is to treat them as limiting cases of extended sources, with the distribution 
tending towards an impulse or delta-function. If this is to be expressed in a coordinate 
system of polar and azimuth angles, one has to take into account the non-uniform 
spacing of coordinates. Consider a collimated source which produces an irradiance E 0 
on a surface oriented orthogonally to the direction ( 6 0 , 0 O ) of its rays. Clearly the 
radiance Lj of this source should be zero except for this direction. The product of 
Dirac delta functions, 5 (8,- — 6q) • 5 (0 ; - — 0q), will be a useful ingredient of the 
formula expressing Lj as a function of the angles. One must insure, however, that th<». 
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irradiance on a surface lying orthogonal to the rays equals Eq, that is, 

n ir/2 

Li sin 0/ ddj dty. 

Clearly this can be accomplished if 

^ = E 0 6(0/ - 6q) 6(0/ - 0q) / sin 0 O 

This is called the “double-delta” representation of source radiance for a collimated 
source. It can also be written in an alternate form using the identity 

6 [/(*) ~ /M - *(* ~ *o) / f'( x o)- 

where f(x 0 ) is the derivative of f(x) evaluated at x — Xq, provided that fix) has an 
inverse in the domain of interest. (This identity can be confirmed by integrating each 
side with respect to x over the relevant interval.) Then, 

Lj = Eq 6 (cos 0/ — cos 0 O ) 6(0/ — 0 O ). 


§8. PERFECTLY SPECULAR REFLECTANCE 

A perfectly specular or “mirror-like” surface reflects light rays in such a way that 
the exitant angle 0 r equals the incident angle 0/ and that the incident and reflected ray 
lie in a plane containing the surface normal. The reflected radiance of a surface patch 
in the direction (0 r ,0 r ) is simply the source radiance in the corresponding reflected 
direction. That is, 

£ r (0 r ,0 r ) = Z./(0 r ,0 r +»r) 

The surface thus forms a virtual image of the source. From the definition of the 
BR.DF, we see that 

L r — J* f r dEj = Jq' f r Lj dQj 

That is, 

n ir/2 

f r Lj cos 0/ sin 0/ ddj dty 

« 

We can satisfy the conditions stated above if we let 

fr.is ~ ^ (0/ — 0 r + *r) / (sin 0/ cos 0/) 
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This is called the “double-delta” form of the BRDF for perfectly specular reflectance 
[43? ^ ide " t,ty me " ti0ned ' n th£ 1351 ■“** ^ “ WHte ^ ^ f0r ” 


fr,is — 2 6 (sin 2 0 r — sin 2 #,-) 6 (0 r — 0/ -4- tt) 


§9. ANALYSIS OF IMAGE-FORMING SYSTEM 


We will now analyze a simple image-forming system (Fig. 6). We assume that the 

device is properly focused; that is, those rays originating from a particular point on the 

object which pass through the entrance aperture are deflected to meet at a single point 

in the image plane. Similarly, rays originating in the infinitesimal area dAg on the 

object s surface are projected into some area dAp in the image plane and no rays from 

other portions of the object’s surface will reach this area of the image. Further, we 

assume that there is no “vignetting,” that is, the entrance aperture is a constant circle 

of diameter d and does not become occluded for directions which make a larger angle 

with the optical axis. The effect of vignetting on image irradiance will be considered 
later. 


The exposure of film in a camera is proportional to image irradiance, E p , and gray- 
levels in a digital imaging system are quantized measurements of image irradiance. In 
order to calculate image irradiance we must first determine the flux d$ L passing 
through the entrance aperture arriving from the patch of area dAg on the object. 


L 


where Q r is the projected solid angle subtended by the aperture. We will also need to 

know the area dA p of the image of the patch, since image irradiance E p is the flux per 
unit area: p 


E p = d$i/dA p 

If 0/ is the angle between the normal on the surface and the line to the entrance 

aperture nodal point, while a is the angle between this line and the optical axis, then, 
by equating solid angles, 


(dAg cos e/)/f 0 2 = (dA p cos a)// p 2 
E p — (fg/fp9- cos a f L r (cos 6 r /cos #/) dw r 

w r 


Consequently, 
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FIGURE 6 : A simple image forming system. Light collected by the lens from 
the surface patch of area dA§ is projected into the image patch of area dAp . 


Here the integral is over the solid angle occupied by the entrance aperture as seen from 
the patch on the surface. If we assume that the lens is small relative to its distance 
from the object, 6 r is approximately the same as 0/, and the ratio of their cosines is 
unity. Furthermore, the reflected radiance L r will then tend to be constant and can be 
removed from inside the integral. The solid angle occupied by the lens as seen from the 
surface patch is approximately equal to the foreshortened area (rr/4) cP- cos a, divided 
by the distance (f q/ cos a) squared. Finally, one obtains the well-known result, 

E p = L r (tr/4) (d/f p ) 2 cos 4 a 

That is, image irradiance is proportional to scene radiance. The factor of 
proportionality is n divided by four times the effective f-number (f p /d) squared, times 
the fourth power of the cosine of the off-axis angle, a. Thus the “sensitivity” of such an 
imaging system is not uniform over an image, but is constant for a particular point in 
the image. Vignetting introduces an additional variation with image position. Ideally, 
an imaging device should be calibrated so that this variation in sensitivity as a function 
of a can be removed. 









Horn &. Sjoberg 



October, 1978 


Other kinds of imaging systems, such as microscopes or mechanical scanners, lead to 

somewhat different expressions. Generally, however, image irradiance is proportional to 

scene radiance in such systems too. At this point we should remember that scene 

radiance depends on properties of the surface layer (BRDF) and the distribution of 
light-sources (source radiance) since 

— Jq . f r Li dQ-i 


§10. VIEWER-ORIENTED COORDINATE SYSTEM. 

So far we have considered directions from the object to the image-forming system 
and to light sources in terms of a local coordinate system with one axis lined up with 
the surface normal. Such coordinate systems will vary in orientation from place to 
place and are thus inconvenient for the specification of global distributions such as that 
of source radiance. A coordinate system fixed in space will be more suitable, 
particulary if one of the axes is lined up with the optical axis (Fig. 7). In this viewer- 
oriented coordinate system we introduce polar angle $ measured from the z-axis and 
azimuth angle 0 measured from the x-axis in the plane perpendicular to the z-axis. 
Here, the z-axis is parallel to the optical axis. Directions to sources of light can be 
given using these two angles. If the sources are far away in comparison to the size of 
the objects being imaged, then source radiance will be a fixed function of these angles 
independent of the point on the surface being considered. 
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§11. the surface normal 

In the local coordinate system the surface normal lies along one of the axes, or 
equivalently, it is the direction corresponding to zero polar angle. In the viewer-oriented 
coordinate system the surface normal will correspond to some direction, say (). 
The corresponding unit vector is, 

n = (cos <f> n sin 6 „, sin 0„ sin $ n , cos $ n ) 

The surface of the object may be specified by giving “elevation” z as a function of the 
coordinates x and y. We can give an expression for the surface normal in terms of the 
first partial derivatives of z with respect to x and y, if these exist. Let the first partial 
derivatives be called p and q. Then the vectors (1,0, p) and (0,1 ,q) are tangent to the 
surface, as can be seen by considering infinitesimal steps in the x and y direction. The 

surface normal is perpendicular to all vectors in the tangent plane and so is parallel to 
the cross-product of these two: 

(It 0, p) X (0,1, q) = (—p, —q, 1) 

. Thus the unit normal can be written 

n = (— p, — q, 1) / yj 1 -f- p 2 + q 2 

* * 

The following results are obtained by equating terms in the two expressions for the 
surface normal: 

sin 6 n = yjp 2 + q 2 / ^1 + p 2 + q 2 
cos 6 n = 1/ yj 1 p 2 •{- ~qj 2 

sin 0„ = -q / yjp 2 + q 2 
cos 0„ = — p / yjp 2 + q 2 

Conversely, 


p = —cos 0„ tan 8„ 
q — —sin 0„ tan d„ 
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§12. RELATIONSHIP BETWEEN LOCAL AND VIEWER-ORIENTED 
COORDINATE SYSTEMS 

In order to calculate the scene radiance, we will integrate the product of the BRDF 

and the source radiance over all incident directions. Since the BRDF is specified in 

terms of the local coordinate system, while the distribution of source radiance is likely 

to be given in the viewer-oriented coordinate system, it will be necessary to convert 

between the two. Given the direction of the surface normal, (6„ ,0„ ), and the direction 

to a portion of the source (8 s ,<p s ), both specified in the viewer-oriented system (Fig. 8), 

we have to find the incident direction (6,-, 0,) and the exitant direction (0 r , 0 r ) both 

specified in the local system. Alternatively, given the surface normal and the incident 

direction we may have to Find the direction to the source and the exitant direction. 

Note that 6 r = 6„ since the exitant ray lies along the z-axis in the direction towards 

the viewer. Further, since we have excluded anisotropic surfaces, we are only interested 

in the difference between 0 r and 0,-. From the relevant spherical triangle (Fig. 9) we 
obtain 



viewer-oriented coordinate system. 
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FIGURE 9: Spherical triangle extracted from previous figure 
and used in derivation of transformation equations between the 
surface-normal, local coordinate system, and the viewer-oriented, 
global coordinate system. 


Cosine formula: 

cos 6j = cos 0, cos 0„ + sin 0, sin 0„ cos (0,-0,,) 

Sine formula: 

sin Qj sin (0 r — 0/) = sin 0, sin (0, —<t> n ) 

Analogue formula: 

sin 0,- cos (0 r — 0/) = cos 0, sin 0„ — sin 0, cos 0„ cos (0, — 0„) 

The Jacobian of the transformation from ( 05 , 0 ,) to (0/,0/) equals 

( 00 , 700 ,) (00,750,) - (00,730,) (00,/30,) = (sin 0,/sin 0,). 

(The Jacobian will be required below when converting a double integral with respect to 
one set of coordinates to one in terms of the other.) The above formulae allow us to 
find the incident direction from the source direction. Quite symmetrically, we can also 
obtain the source direction from the incident direction: 


Cosine formula: 

cos 0, = cos 0,- cos 0 r + sin 0,- sin 8 r cos (0 r — 0,-) 

Sine formula: 

sin 0, sin (0, — 0„) = sin 0,- sin (0 r — 0,-) 

Analogue formula: 

sin 0, cos (0, — 0„) = cos 0,- sin 0 r — sin 0,- cos 0 r cos (0 r — 0,*) 
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The Jacobian of the transformation from (0,-, 0,) to (8 S ,<t > s ) equals 

(d6 s /ddj) (00, fb<t>i) — (d<t) s /d6j) (dd s /d<t>,) = (sin 0,/sin 8 S ). 

§13. SCENE RADIANCE 


It follows from the definition of the BRDF that reflected radiance can be written as 
the integral 

L r = Jq' f r Lj dSlj = j*^' f r Lj cos 8} du>j 

$' 

Using polar and azimuthal angles this becomes 

n ff/2 

fr(8j,<t>j", 8 r ,0^.) Lj(8 s ,<t> s ) cos 6j sin 8j d8j d<t>i 


Here we integrate over all possible incident directions (0,-,0/) and calculate source 
direction (0,,0 5 ) from the given surface normal (0„,0„) and the incident direction. 
The inner integral has the limits 0 to 7r/2 for 0,-, corresponding to directions within the 
hemisphere visible from the surface. The integration can be extended to the full sphere 
of directions if the integrand is forced to be zero when 8j lies between ir/2 and n. This 
can be accomplished by replacing cos 0/ by max [0, cos 0/ ]. Hence 

f r Lj max[0, cos 0/] sin 0 ; - d8 t dfy 


Since the integral now is over the full sphere of directions, it can be easily rewritten 

using any other set of polar and azimuth angles. Using the viewer-oriented coordinate 
system, for example, we obtain 


L r = X rJo fr L ‘ max f°’ 008 6 A sin e s d6 s d<t>s 

That is, 


L r (8 n ,<t> n ) — frWi’ty'* Or’Qr) Lj(8 s ,<t> s ) max[ 0, cos 0,] sin 8 S d8 s d<t> s 

Here we integrate over all possible source directions (0 S ,<£> S ) and calculate incident 
directions (0/,0/) from the given surface normal (0„,0 n ) and the source direction. We 
now have two convenient forms for the calculation of scene radiance. We proceed to 
calculate reflectance maps for a few simple combinations of BRDF and distributions of 
source radiance. 
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§14a. COLLIMATED SOURCE, LAMBERTIAN REFLECTANCE 
For a lambertian reflector, f r = 1/jt. For a collimated source, 

L { tm Eq 6(0, — 0 O ) 6(0, — 0 O )/sin 0 O 

where Eq is the irradiance measured perpendicular to the beam of light arriving from 
source direction (0 o ,0q)- Substituting into the second form of the expression for scene 
radiance above, we obtain 

n ir 

(E 0 /n) 6(0, — 6 0 ) 6(0, — 0 O ) max[0, cos 0,] (sin 0,/sin 0 O ) d$ s d<f> s 
This is equal to 

# 

L r = {Eq/ tt) max [0, cos 0/ ] 
where 0/ — 0,- when 0, = 0 O and 0, = 0 O . In this case, 

cos 0/ = cos 0 r cos 0 O + sin 0 r sin 0 O cos (0 O — 0„) 

and 

cos (0 O — 0 fl ) = cos 0 O cos 0„ — sin 0 O sin 0„ 

To obtain the reflectance map, scene radiance as a function of surface gradient, we can 
substitute expressions in p and q for these trigonometric expressions. The result is 


R(p , q) = (E 0 / w) max 


0 , 




yj l+p 2 +q 2 yjl+pQ 2 +qQ 2 


where 


Pq = —cos 0 O tan 0 O 
go = —sin 0o tan 0 O 

A surface with gradient (p 0 , q 0 ) is normal to the direction of the incident light rays. 
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§14b. UNIFORM SOURCE, LAMBERTIAN REFLECTANCE 

A uniform source has constant incident radiance. Let Lj = Lq. Again, for a 

lambertian reflector, f r — 1/tt. Substituting into the first form of the expression for 
scene radiance, we obtain 



‘it Orr/2 

nJo (W*) 005 fy sin Mi d<t>i 


This becomes 

r*P- 

L r = L 0 Jo Sin 26 i d6 i = A) 

Not surprisingly, the reflected radiance is independent of the surface orientation in this 
case. 


§14c. HEMISPHERICAL UNIFORM SOURCE, LAMBERTIAN 
REFLECTANCE 

A hemispherical uniform source is described by 

£/(0* «0$) = A) f° r 9 $ < n /2 

£/(^i0j) = 0 for $ $ > ir/2 

In order to evaluate the double integral for scene radiance, it is necessary to know the 
value 0/ of the incident angle 0,-, which corresponds to the “horizon” of the “sky,” that 
is, 6 S — rr/2. From the coordinate transformation equations one can easily see that 


cot 0/ = —tan 0 r cos (0 r — 0,) 

For <t> r — n/2 < 0,- < 0 r + *r/2, the horizon cutoff only occurs for 0/ > ir/2 and 

can be ignored, since the inner integral is from 0, = 0 to 0,- s= »r/2 only. For the 
other half of the range of 0 /f this cutoff must be considered. Now, 

n n/2 

/ r Z,-cos 0,-sin 0,-</0,-</0,- 



r* rmin[0/,7r/2] 

L r — (Lq/ *) J_ w J q cos 0/ sin 0,- d6,- dty 


Letting 0 = 0/ — 0 r and 0' = 0 ; - — 0 r + tt, we can split the outer integral. 
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pr/2 rn/2 pr/2 p$( 

L r — (A)/”') J_ n/2 J 0 cos 0, sin0 t - ddj d<f> + (Lq/ n) j ^ J q cos 0, sin 0, dOj d<f> 


r*n 

Jo c 


cos 0j sin 0j d0j = 1/2 


so the first term is simply Lq/ 2. Next, note that 

f 0! 

J 0 cos 0j sindj ddj d0j = (1 — cos 20/)/4 
where, since cot 0( = —tan 0 r cos (0 r — 0 ; ), 


sin 2 0-/2 


sin 2 = 1/[1 -f- tan 2 0 r cos 2 (0 r — 
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FIGURE 10: Cross-section through uniform hemispherical source and surface 

element, illustrating horizon cutoff and portion of extended source not visible 
from surface. 
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Z~\ 


The second term thus becomes 


r /2 

(Lq/Itt) J 1/[1 + tan 2 6 r cos 2 <f/] d<$> 


which equals 


(L 0 /2tt) [cos 6 r tan-1 (cos 6 r tan <t>)]^ = (Lq/2) cos 8 r 


Adding up the two terms we finally get 


Lr^n'^n) ~ (Lq/2) 0 + cos 0 n ) = Lq cos 2 (0 fl /2) 


This is the result found by Brooks [44]. From it the reflectance map can be found 
immediately: 


R(P > 7) — (Eq/2) (1 -f- 1, 


<7 2 ) 


§15a. COLLIMATED SOURCE, SPECULAR REFLECTANCE 
For specular surfaces, 


fr — *>(€>,■ — 6 r ) 6(0/ — 0 r + 7 r) / (sin 0/ cos 0,) 

Using the source radiance from §14a above, and the first form of the expression for 
scene radiance, we obtain 

/**■ /*jr/2 


a 


That is, 


where 0 / and 0 / 
0 / — 0 r + rr. Usii 
0 / = 2 0 r and 0 / = 


L r = Lq 6 ( 0 / — 0 O ) 5 ( 0 / — 0 O ) / sin 0 O 

are the values of 8 S and 0 5 corresponding to 0 / = 0 r and 

5 the equations for the coordinate transformations, one finds that 
0 rt . Thus, 

L r = Lq 6 ( 20 r — 0 O ) 5 ( 0 „ — 0 O ) / sin 0 O 


and finally, 


L r (8 n , 0„) — (Lq/2) 6(0„ — 0q/2) 5(0„ — 0q) / sin 6t 
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To express this as a function of p and q we have to remember that 

6 £/(■* . y) — f(x o, j> 0 )] 6 [g(x, y) - g(x 0 , y 0 )] = 6 (x — Xq) 8(y — y 0 )/J(x Q , y 0 ) 

where J(x, y) is the Jacobian of the transformation from (x,y) to (f, g), (provided this 
transformation has an inverse in the region of interest): 

J( x - y) — (df/dx) (dg/By) — (df/dy) (dg/dx) 

The Jacobian of the transformation from (p,q) to (B„, <f> n ) is 

i 

— l/[ 'Jp 2 + q 2 (1 + p 2 + q 2 ) ] 

Let 


Pi — —cos 0 O tan 6q/ 2 
q\ — —sin 0 O tan B 0 /2. 

Then, noting that sin Bq ~ 2 sin (Bq/2) cos (6q/2), one can write 

sin Bq — 2 yjp\ 2 + q\ 2 / (1 4 - p\ 2 + q\ 2 ) 

and therefore, 

R(p,q) = (E 0 /4)6(p - p { ) 6(q - q { ) (1 + p { 2 + ^,2)2 


Thus, a surface element with gradient (/>j»<7j) is oriented to specularly reflect the 
collimated source towards the viewer. This gradient can be related to the gradient (jfn, 
q$) introduced earlier: 

Pi ~ PO (V 1 + Po 2 + Bo 2 — 0 / (Pq 2 + qo 2 ) 

Q\ ~ % (V 1 + Po 2 + Bo 2 “ 1) / (Po 2 + ^o 2 ) 

When the point (Pq.^q) is not far from the origin, then (p\,q{) is approximately 
midway between the origin and (p 0 , q§). 
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§15b. UNIFORM SOURCE, SPECULAR REFLECTANCE 

It is easy to see that for a specular surface under a uniform source, the scene 
radiance will be constant and equal to the source radiance. 

L r — Lq 

This is the same result as the one we obtained for the uniform source and lambertian 

reflectance. Thus a diffuse surface appears just as bright as a specular surface if both 

are viewed with uniform illumination. In fact, all surfaces reflecting the same fraction, 

p say, of the total incident light will appear equally bright under this illumination 
condition. 


§15c. HEMISPHERICAL UNIFORM SOURCE, SPECULAR REFLECTANCE 

In this case, 

= Lq for < it/A 
L r (O n ,<f> n ) = 0 for 0 n > it/A 

The reflectance map is 


R(P ,q) = Lq for p 2 + ql < i 

R(p ,.q) = 0 for p 2 -f- <p- > 1 


SUMMARY AND CONCLUSIONS 

We have shown that image irradiance is proportional to scene radiance and that 

scene radiance depends on surface orientation. The reflectance map gives scene 

radiance as a function of the gradient. It can be calculated from the bidirectional 

reflectance-distribution function (BRDF) and the distribution of source radiance. 

Several special cases were worked out in detail. Each could have been developed more 

easily by a direct method, but was obtained from the general expression for scene 

radiance to illustrate the technique. The general expression allows one to find the 

reflectance map even if the source radiance distribution or the BRDF is only given 
numerically. 
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